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Abstract: The note presents an algorithm for the Smarandache's function computation. The 
complexity of algorithm is studied using the main properties of function. An interesting inequality is 


found giving the complexity of the function on the set {1,2,...,n}. 
1. Introduction 


In this section, the main properties of function are reviewed. The Smarandache’s 

function notion reported for the first time in [1]. The main results concerning the 
function can be found in [1], [2]. 
The function S:N—WN defined by S(n)=min{k{k!in} is called Smarandache’s 
function. This concept is connected with the prime number concept, because using the 
prime numbers an expresion for the function is given. The important properties that 
are used in this paper, are showed bellow. 


1. For all neéN, the inequality S(n) <mis true. When n is a prime number, the 


equality is obtained . 
2.1f mn=ph-ps..-p™ is the prime number decomposition then 
S(n) = max{S(pj"),S(py"), (pe )} (1) 
_ 3. The inequality S(p*) < p-k is true, if p is a prime number. (2) 


A lot of conjectures or open problems related ot the Smarandache’s function appear in 
the number theory. Several such problems have been studied using computersand 


reported in relevant literature, e.g.{3], [4]. Using the computer 


2. An Algorithm for the Smarandache function 

In the following, an algorithms for computing the function S is presented. The main 
idea of the algorithm is to avoid the calculations of factorial, because the values of 
n!, forn>JO are avery big number and cannot be calculated using a computer. 

Let (x, ),., @ sequence of integer numbers defined by x, = k! mod n,(V)n>0. Using 
the definition of sequence , the following equations can be written: 

° x, =1 (3) 
© xX,,, =(k +1)! mod n=(k+1)k!modn=(k+1)x, modn. (4) 
Based on the linear equation (4), S can be calculated § as 
S(n) = min{k|k!in} = min{k|x, = 0}. 

The algorithm for S(n) performs the computation of X,,X,...,X, until the 0 value is 


found. The PASCAL description of this algorithm is given bellow. 


function S(n:integer):integer; 
var 
k,x:integer; 
begin 
x:=1;k:=1; 
while x0 do 
begin 
x:=x*(k+1) mod n; 
k:=k+1; 


end; 


An analysis for the complexity of algorithm is presented in the following. The work- 


case complexity and the average complexity are studied. 
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Theorem 2.1 

If n= pi'-p;*-..-pe" then the complexity of the algorithm for computing S(n) is 
O(max{p,k,, D2 ---sPmkm})- 

Proof 

The algorithm computes the value S(n) generating the sequence X1yXq5-++sXgiqy- The 
number of operation has the complexity O(/S(n)). 

Based on (1) and (2), the following inequality holds 

S(n) = max{S(p;"),S(p3"),---.S(par)} S$ max {py ky, Pyky 1 Paka} - (5) 
Therefore, it can “be concluded that the complexity of computing S(n) is 


O(max {p,k, »P,k, cs 2 


Other aspect of complexity is given by the average operations number. Assume the 
value S(k) is generated, where k is a number between | and n. This value can be 
computed with S(k) operations. Therefore the process takes S(1) operations for the 
value S(1), S(2) operations for the value S(2), ...,.a.s.0. The average of the numbers 


operations is 5=-Y Si and gives an other estimation for the complexity of 


ist 


algorithm. 

In the following, a possible upper bounds for S are sought. Obviously, S verifies the 

simple inequality 

=_1< be et A 

S=-—) S@<—-—) i=—=—n+-. 6 
n 2 @) gee eae ae (8) 


Inequality (6) can be improved using the strong inequalities for S/i). 


Lemma 2.1. 


If >2 is an integer number the following inequalities hold 

1. S(2i)+S(2i+ Dit]. (7) 
2. S(2i-1)+S(2i)Bi-1. (8) 

Proof 


SO 


Assuming i>2 folow S(2i) <i. Applying this result both inequality are true. 
Based on lemma 2.1, we can found an upper bound for the S better than in (6). 


Theorem 2.2 


If n> J is a integer number then the inequality S= “> sts < on + zt + at is true. 


i=] 


Proof 
Two cases are considered to prove the inequality. 


Case 1: n=2m 


nS = y S(i) = 3 [S(2i — 1) + S(2i)] = S(1) + S(2) + S(3) + S(4) + y [S(2i-1) + S(2i)] 


ix} i=l i=3 


Based on (8) it can be derived that 


nS = 945° [s(2i—1) + 5(2] <9 +5 Gi- l= 245° Gi-)= 2452 G+)-2- 


ix} iz} i=} 


spate ayy BO ee pa: 
29-2 i aa 
Dividing by 7, we obtain the inequality S<eneies, (9) 
n 


Case 2: n=2m+1 


nS = 3 S() =S() + Sines + S(2i +1)] = S(2) + S(3) + S(4) + S(5) + “S*ts(2 + S(2i+1)] 


i=l i=} 


Using (7), it is found 


f (a-2 (a-D/2 (a-)2 caren ewe n-l 
nS=14+ S'[S(2i)+S(2i+)] S14+ VGit)=3+ Yi Git) =3+= oe 
ia} i=} i=l 
= 3+ =(n? -1) Here Re va Thus, S22 7el eg! (10) 
2 2 8n 
From (9) and (10), itis found S <—n + min{ = a4 isons rote 
n 


3. Final Remarks 


$1 


1. Based on theorem 2.2 we can say that the average operations number for 


3 2 
computing the Smarandache’s function is less than =” + : +=. 
n 


2. The upper bound =n + - +2 improves the previous bound zn + : : 

3. The improving process can be extended using other sort of inequalities give the 
prime numbers 2 and 3. A lemma as similar as lemma 2.1 finds the upper bounds 
the sum of sixth consecutive terms of Smarandache’s function. 

4. Using the algorithm for computing the function S, the Smarandache’s function can 
be tabulated. The values S(1),...,S(n) for all n <5000 can be found. The algorithm 


should be reviewed to be able to compute the Smarandache function for the big 


numbers. 
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